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Abstract 

Quantum operations are usually defined as completely positive (CP), trace preserving 
(TP) maps on quantum states, and can be represented by operator-sum or Kraus repre- 
sentations. In this paper, we calculate operator-sum representation and master equation 
of an exactly solvable dynamic of one-qubit open system in layered environment . On 
the other hand, we obtain exact Nakajima-Zwanzig (NZ) and time-convolutionless (TCL) 
master equation from the maps. Finally, we study a simple example to consider the re- 
lation between CP maps and initial quantum correlation and show that vanishing initial 
quantum correlation is not necessary for CP maps. 

Keywords : Kraus representations, Quantum Correlation, Nakajima-Zwanzing 
and Time-convolutionless Master Equation 

1 Introduction 

In quantum information theory, it is widely accepted that any physical process can be described 
by a quantum operation or quantum channel which is CP (this means that not only map $ : 
Bi — > B 2 on c*-algebra is positive, but also the combined operation <&®In : Bi <g> N — > B 2 <8> N 
for all dimensions N is positivity preserving ) and TP maps (this means that Tr$(p) = Tip for 
all trace class operator p) between spaces of operators pQ- [5]. These maps play an important 
role in the description of nonunitary time evolutions of open quantum systems that interact 
with an environment. On the other hand, the environment is classified as Markovian with no 
memory effect and non-Markovian with memory effect [6]-|19j. 

It is well-known that for closed quantum system, its time evolution can be described by a 
unitary operator. However, for an open system, the time evolution is not necessarily unitary. 
The evolution of an open system is usually described by the Kraus representation [20J. The 
Kraus representation of an open system is usually constructed by considering a large closed 
system. Let's assume an interaction between the open system denoted as (s) and the envi- 



ronment (b). This environment is a quantum system with the Hilbert space of an arbitrary 
dimension. The whole (s) + (b) system evolves unitarily. In most of the studies on dynamics 
of open systems, it is assumed that the open system and its environment are at the initial 
moment of their joint evolution factorized|21[ 122], that is they are described by the density 
operator of the form 

Psb = p s ® Pb, 

where p s is the initial state of system and p b is the initial state of the environment. As the 
combined system is a closed one, its evolution is unitary, 

p sb (t) = U sb p sb {0)Ul, (1.1) 

where U sb is the unitary operator. The interested system, as an open one, than evolves in the 
following way 

Ps(t) = tr b {U sbPsb Ul} = tr b {p sb (t)}. (1.2) 

The above equation is obtained by doing partial-trace on bath. Krausj2ni E3] and Choi [21] 
showed that if the above equation can be equivalently expressed in the form 

p s (t) = $(p s (0)) = ^M^(t)p s (0)Mt,(t), (1.3) 

where the map $ is described by $ : p s (0) — > p s (t) and M^{t) satisfy the following equation 

Y,Ml(t)M, u (t) = I 1 (1.4) 

/«/ 

it is said that evolution of p s (t) has the form of the Kraus representation and the above equation 
is condition of CP map[25j. 

In this paper, we consider a model of a spin star configuration which consists of N+l 
localized spin-| particles. One of the spins is located at the center of the star, while the 
other spins on concentric circles with different radii surrounding the central spin layer by 
layer |26j. The difference in radius is because the coupling coefficients between layers spins and 



the central spin is taken differently. At first, we obtain Kraus representations and then arrive 
master equation with the use of Kraus representations and show that they satisfy Eq.(4). 

Maps are one of the main bases to obtain Kraus representations. In addition, they are 
effective in obtaining exact NZ and TCL master equation[27|. NZ and TCL are used when 
the exact solution of the master equation is a difficult task, therefore we can approximate 
the master equation to solve [281 |Sj. These techniques were introduced by Nakajima (1958) 
and Zwanzig (1960) and independently by the Brussels school (Prigogine, 1962) [8j. They are 
widely used in non-equilibrium statistical mechanics and have played an important role in 
various fields of physics (e.g. in the theory of transport phenomena and relaxation processes, 
in quantum theory damping, in super-radiance and laser theory) [29] . Details of the application 
of this method in the various fields and numerous references are given in books and reviews 
articles of H aken [301 EI], Agarwal [32] and Haake [33]. By considering the above model and 
exploiting the knowledge of the exact unitary evolution, and the reduced dynamics, as well as 
a suitable matrix representation of the dynamical maps, we can exhibit the exact TCL and 
NZ equation of motion. 

The paper is organized as follows. In Sec. (II), we will explain about how to obtain map and 
Kraus representation and master equation. Exact NZ and TCL master equation are discussed 
in Sec. (III). Finally, in Sec. (IV), we consider a simple example that shows vanishing quantum 
discord is not necessary for CP maps. 

2 Map, Kraus representation and Master equation 

2.1 Background 

The first step in relating a general, non-Markovian master equation p(t) = A[p(t)] to a cor- 
responding CP map p{t) = $t[p(0)] relies on expressing this relationship in matrix form. We 
need to write the matrix form to the basis of {H 7 ^}, that the basis denotes every convenient 



orthonormal basis set for the Hermitian operators on the Hilbert space, i.e. 

wl = W a , 

tr[W a W b ] = 5 ab . 

Now, according to the introduced basis, for any maps $ (though we will be concerned mainly 
with CP maps below) we can express both p and $(p) in matrix form as follows [34|, 

$(p) = (Fr) T W, (2.5) 

where the evolution matrix F is the following 

F kl :=tr[W k $(Wi)], (2.6) 

and also vector r is defined as 

n :=tr[Wip\. 

Indicating time-dependence of $ and F via $t and F(t), and defining 

p(t):=$ t [p(0)]. (2.7) 

Time evolution of p can be expressed as 

p = [F(t)r(0)} T W. 

Now, suppose that p(t) satisfies a master equation of the form 

p = A(p). (2.8) 

Note that the linear map A may be time- dependent. Now, this master equation can be rewrit- 
ten as 

p = [L(t)r(t)] T W. 



Note that L (and F) are necessarily real, but not symmetric in general. Comparing the results 
of the last two paragraphs gives 

Fr(Q) = Lr{t) = LFr(0). (2.9) 

Hence, since by linearity this equation must hold for all vectors r(0) (whether or not they 
correspond to a density operator). 

Choi has shown that complete positivity of $ is equivalent to the positivity of the par- 
ticular matrix S [351 EH], and that the Kraus decompositions of $ are related to the outer 
product decomposition of S. The master equation gives us matrix L, which in turn gives us a 
matrix F, characterising the linear evolution map $f In order to proceed from F to a Kraus 
decomposition, all we need to do is finding S from F, and then diagonalise it. In fact, our 
construction of S from L allows us to determine whether or not, for a given proposed master 
equation, the corresponding map is completely positive. Choi further demonstrated that $ is 
completely positive if and only if S is positive, i.e., S > [351 ES]. The Hermitian property 
implies that one can always decompose S as a sum of outer products, i.e., 

S = J2v t (t)vHt) (2.10) 

i 

for some set vectors{F(z)} that we take {V(i)} as follows 



v(i) = V[x~\n(i), (2.ii) 

where II coefficients are eigenvectors of matrix S and also the Kraus operators are given by 

M< := £V(i) a &. (2.12) 

a 

Noting that the j3 a form an orthonormal (non- Hermitian) basis for the operators on the Hilbert 
space [51] . 



2.2 Map, Master Equation and Kraus representation for spin star 
model for one-qubit system with layered environment in a chan- 
nel 

We obtain all the issues described in the previous section for the model which at first was 
mentioned. The Hamiltonian of the XX central spin model that is composed by a localized 
spin, hereafter called central spin, coupled to N spins of bath layers with the different coupling 
constant a, takes the form [26] 

# = 2(a+H_+a_H + ), (2.13) 



where E± is denoted as 



/ , a »J+i 



n 



(2.14) 
(2.15) 



where a± and J± are the Pauli operators referring to the central spin and the surrounding 
spins respectively , and J± is denoted as 



So, evolution density matrix of system is as follows 

/ 



(2.16) 



Ps(t) 



\ 



l + u 3 (0)/ 3 i>-(0)/ 12 
u+(0)/i 2 l-v 3 (0)/ 3 



/ 



(2.17) 



where we have introduced the functions|26| 



f 12 (t) = tr B {cos[2thi(ai, . . . , a n )] cos[2t/i 2 (ai, . . . , a n )\ ® 2 I B } 



(2.18) 



and 



f 3 {t) = tr B {cos[2t/»i(ai, . . . , a n )] ® 2"^} 



(2.19) 



where hi(ai, a 2 , • • • , <x n ) and h 2 (a\, a 2 , ■ ■ ■ , a n ) are 

/ii = hi(ai,a 2 , 



i C*r. 






/i 2 = /i 2 (ai,a2,---,ar. 



\ 



£«W 



(2.20) 



(2.21) 



To obtain the master equation, the matrix F firs must be determined, and checked as whether 
it is invertible. The basis {W a } is chosen as W a = -j^Pa, where a a are the Pauli operators. 
One finds 



$(Wq 



$(Wi 



v/2 



1 



1 



\ 



/: 



12 



V2 



/l2 



* (Wa) = -L 


-i/i 5 
v i/12 


/ 


$(w 3 ) = 4 

V 


2 1 o A , 





The matrix F follows via Eqs.(5) and (6) as 

/ 



F 



10 

/12 o 

o / 12 

V o / 3 J 

The solution for L follows via Eq.(9)as 

/ 



L = FF 



-i 




0^00 

oo^o 



/12 



/£ 



\° ° ° !/ 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



(2.27) 
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Now, to obtain the master equation in the form p = A t (p), we calculate the Choi matrix R in 
Ref.pg], 



R 



2/3 


f{2 

h-2 








\ 


f{ 2 
hi 


2/3 












i n 



n 


fk 

2/3 

n 




fk 





(2.28) 



The master equation immediately follows via equation as 



A t (p):=J2Rab(t)/3 a p(3l 



Pit) = At[p(t)} 



a I) 



J 3 , J 12 u 1 J 3 / i 



(2.29) 



'2/ 3 /l2^ " ^ " " J 2/ 3 

Finally, we obtain Kraus representations via using of Eqs.(ll) and (12). At first, we calculate 
matrix S as follows 



/ 



s = i 

2 



l + 2/i 2 + / 3 

l-/ 3 

OI-/3 



\ 



V 











(2.30) 



1 - 2/ 12 + h ) 



Now, we obtain Kraus representations 



U = v /2(l + 2/ 12 + / 3 )Wb, 



M 2 = A / -(1-/3)^1, 



-U 5 V^(i-/ 3 )w 2 , 



M 4 = M /-(l-2/ 12 + / 3 )W 3 , 



(2.31) 

(2.32) 
(2.33) 
(2.34) 



that they satisfy Eq.(4). 
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3 exact Nakajima-Zwanzig and time-convolutionless mas- 
ter equation 

3.1 Background 

With the aid of the knowledge of the exact time evolution, and using the representation of 
maps in terms of matrices, we now explicitly obtain two kinds of exact equations of motion for 
the reduced system's dynamics. We first consider a master equation in differential form with 
a generator local in time, that is, the TCL master equation. Assuming the existence of such a 
generator K TC L(t), it should obey the equation [2] 



p{t) = K TCL (t)p(t), (3.35) 

which, due to p{t) = <3>(£)p(0), is satisfied upon identifying 

K TCL = ^^(i) (3.36) 

or, in terms of matrices, 

Ptcl = Ht)F-\t), (3.37) 

this expression holds when F is invertible. 

Now, we can obtain the NZ master equation according to the TCL master equation when 
memory kernel exists. In this case the memory kernel Kj^zif) should obey the convolution 
equation 

p(t) = (K NZ op)(t), (3.38) 

so that in view of Eq.(7), one has the relation 

K NZ (u) = ul- $ _1 (V), (3.39) 

where the hat denotes the Laplace transform, and therefore in matrix representation, 

P NZ (u) = uI-F~\u). (3.40) 
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3.2 Nakajima-Zwanzig(A^Z) and Time-Convolutionless(TCL) for spin 
star model for one-qubit system with layered environment in a 
channel 

The technique used in the previous subsection to obtain the TCL and NZ equation of motion 
for a model whose evolution is known, by exploiting the representation of maps in terms of 
matrices, is applicable for a detailed study of Model in Sec.(IIB). 

Now, starting from Eq.(17) and exploiting the same strategy used in Sec.(IIIA), one im- 
mediately obtains, for the matrix representation of the TCL generator, the expression 

/ n n n n \ 



TCL 



(t) 





7i 

7i 

V 72 / 



(3.41) 



where 71 and 72 are 



7i 



-2(/ii tan[2hit] + h 2 tan[2h 2 t}) , j 2 = -4/iitan[4M]. 



And one can also determine the expression of the NZ memory kernel, whose Laplace trans- 
form is given by 



Pnz{u) 



/ \ 

771 

771 

V rj 2 J 



(3.42) 



where rji and rj 2 are 



2 N [u 4 + 8(h 2 + h 2 )u 2 + 16(/i? - hi) 2 } 2 N [u 2 + 16h 2 } 
Vi = u o , . />2 , ,0, ,V2 = U 



u 3 + A{h\ + h\)u u 

So the master equation in operator form reads [27], 

11111 
k T CL{t)p = --r] 2 {(J+P°- - 2^ a - a +-' PS) ~ 2 7l2 ^ a - pa+ ~ 2^ a+a - ,p ^ + 4^ 2 ~ 2r h)( cr 'P <T » ~ P)> 

(3.43) 



12 
and 

K NZ {u)p = --l2{v+pv- ~ 2^- a +' P}) ~ 2^ a -P a + ~ 2^ a+a - 7 p ^ + 4^ 72 ~ ^(vzPVz ~ p)- 

(3.44) 
It immediately appears that, given the time evolution map from the relation Eqs.(41) and (42), 
one can directly obtain the generator of the master equation in TCL form or the memory kernel 
for the NZ form, respectively, without resorting to the evolution of the whole perturbative 
series. Obviously, given the exact time evolution, one dose not need the equations of motion. 

4 Initial quantum correlation and CP maps 

4.1 Background 

Recently, a relation between CP maps and quantum discord has been put forth [37J EH] • quan- 
tum discord, first introduced by Ollivier and Zurek [39], that captures the difference of two 
natural quantum extensions of the classical mutual information, can be used as a measure 
of the quantum correlations. Although quantum discord is equal to the entanglement for 
pure states, it includes the quantum correlation which are contained in mixed states that 
are not entangled. In [37], it is shown that if the initial system-bath state has vanishing 
quantum discord, then the dynamics of system can be described by a CP map. In [38J Sha- 
bani and Lidar made a strong claim: The reduced dynamics of a system is completely-positive, 
for any coupling with the bath, only if the initial system-bath state has vanishing quantum discord 
as measured by the system. 

Here we show that underlying assumption in Ref . [38J limit the generality of the constructed 
reduced dynamics. To do this we consider two different paradigms to ask the question: given 
any unitary evolution for the system-bath, can we define a CP map taking a family of initial 
states {p s (0)} to final states {p s (t)}7 
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4.2 Vanishing quantum discord is not necessary for CP maps 

Now, we prove our claim that vanishing quantum discord is not necessary for CP with notice 
to below example. In this example, first we obtain Kraus representation with map and Eq. (3) 
and then we prove that map is CP with notice to Eq.(4). On the other hand, we compute 
quantum discord and it shows that we can choose conditions for nonzero quantum discord. 

We choose the same model as that Ref.fTB] has used. That is to consider a combined system 
composed of two spin-1/2 subsystems with the interaction Hamiltonian 



H sb = cr x ®-(I-a z ) + I® -(/ + a g ), 



(4.45) 



where a x and a z are Pauli spin operators. In this model, the first qubit plays the role of 
the open system while the second qubit plays the role of the environment. The interaction 
described by the Hamiltonian corresponds to the well-known controlled- NOT gate 
The unitary evolution operator is given by U s f,(t) = e~ iHabt , explicitly 

N 

cost — isint 

e _it 

-isint cost / 



/ e -it 



U ab (t) 





V o 

So, first we obtain p s b(0) as follows 

Psb (0) = -(I s ® I b + J2 *°i ® ° 

and also we can obtain p s (0) with trace of p s b(0) as follows 

1 



(4.46) 



b ) 



(4.47) 



Ps(0) 



-I. 



(4.48) 



where ais are Pauli matrices. From Eqs. (2), (10) and (11), we get the density matrix of the 

system 

l/ 2 + c 3 (l -cos(2t)) -ic 3 sin(2t) 

ic 3 sin(2t) 2 - c 3 (l - cos (2t)) / 



\ 



P,(t) 



(4.49) 
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So, maps can be written, with notice to Eqs (7) and (49), as 

/ 

$(W X 



2v/2 



V 



2 + c 3 (l - cos (2£)) -ic 3 sin (2t) 

ic 3 sm(2t) 2 - c 3 (l - cos (2t)) 

$(Wi) = $(W 2 ) = $(W 3 ) = 0. (4.50) 

Now, we can obtain matrix of F and then matrix of S by having maps. On the other hand, we 
can arrive Kraus representations with have matrix of S as follows 

Mi = ^r(t) sin (t)W + l -T{t) cos (t)W 1 - ^W 3 , (4.51) 

m 2 = -ir(t) cos (t)iy + ^r(t) sin (t)wi + ^w 2 , (4.52) 

M 3 = ir(t) sin (t)W + |r(t) cos (t)W 1 + ^W 3 , (4.53) 

M 4 = ^r(t) cos (t)W - l -T{t) sin {t)W l + ^W 2 , (4.54) 

where T(t) = yl — c 3 sin (t) and Kraus operators {Mj} satisfy Eq.(4) which expresses that 
our choose map is CP. 

However, given the quantum discord calculated in Ref.[40j and Eq.(4), we can conclude 
that according to Kraus representations obtained, Eq.(4) will always be satisfied regardless 
of the values of c 3 and the values of quantum discord and it is clearly violation of claim of 
Ref.j38]. So, we understand that vanishing quantum discordC is not necessary for CP map. 

5 conclusion 

In this paper, we have started from obtained maps of spin star model for one-qubit system with 
layered environment and then by use of these maps, computed the matrix F, With the matrix 
F, first we have obtained the Choi matrix S and then by use of matrix S, have arrived Kraus 
representation and master equation that clearly, we understood that Kraus representation 
satisfy Eq.(4). It expresses that maps are CP. On the other hand, by use of the matrix F we 
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have arrived matrix representations of the TCL and NZ generator and then with use of matrix 
representations we have obtained operator representations of the TCL and NZ generator. In 
the end, with the two subjects expressed (satisfies the equation (4) and the calculated quantum 
discord), we have proven that vanishing quantum discord is not necessary for CP. 
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